Scanning tunneling spectroscopy and Dirac point resonances due to a single Co 

adatom on gated graphene 
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Based on the standard tight-binding model of the graphene 7r-band electronic structure, the 
extended Hiickel model for the adsorbate and graphene carbon atoms, and spin splittings estimated 
from density functional theory (DFT), the Dirac point resonances due to a single cobalt atom 
on graphene are studied. The relaxed geometry of the magnetic adsorbate and the graphene is 
calculated using DFT. The system shows strong spin polarization in the vicinity of the graphene 
Dirac point energy for all values of the gate voltage, due to the spin-splitting of Co 3d orbitals. 
We also model the differential conductance spectra for this system that have been measured in the 
scanning tunneling microscopy (STM) experiments of Brar et al. [Nat. Phys. 7, 43 (2011)]. We 
interpret the experimentally observed behavior of the S-peak in the STM differential conductance 
spectrum as evidence of tunneling between the STM tip and a cobalt-induced Dirac point resonant 
state of the graphene, via a Co 3d orbital. The cobalt ionization state which is determined by the 
energy position of the resonance can be tuned by gate voltage, similar to that seen in the experiment. 



I. INTRODUCTION 

Graphene, a single atomic layer of graphite, has at- 
tracted a great amount of attention because of its fun- 
damental interest and promising applications [1-4]. One 
of the remarkable properties of pristine graphene is the 
linear dispersion relation of its electronic spectrum near 
the Dirac point [5, 6]. It has been demonstrated theoreti- 
cally that chemical adatoms, whether magnetic [7-16] or 
nonmagnetic [16-29], adsorbed on graphene can induce 
resonant states in the vicinity of the Dirac point energy, 
and these states may give rise to strong scattering of elec- 
trons in the graphene. However, to our knowledge, there 
has as yet been no direct experimental observation of 
these Dirac point resonances. The presence of chemical 
adsorbates, atomic or molecular, on graphene strongly 
affects electronic properties and makes graphene-based 
devices suitable for chemical sensing [29-36]. 

Furthermore, the chemical potential of graphene elec- 
trons can be easily tuned by a gate voltage, and thus, the 
adsorbate ionization state can be externally controlled 
[12, 13]. Using single and few-layer graphene films on 
top of a Si02 substrate, it has been demonstrated that 
a graphene field-effect transistor can switch between two 
dimensional electron and hole gases by changing the gate 
voltage [37]. It has been shown that, using an elec- 
trochemically top-gated graphene transistor, even much 
higher electron and hole doping can be induced than by 
standard Si02 back-gating [38]. On the other hand, an 
insulating state with large suppression of conductivity 
can be created in bilaycr graphene field-effect transistors 
by using a top gate in addition to the global back gate 
[39] . When an electric field was applied perpendicular to 
the bilayer, a gap of less than 10 meV at low tempera- 
tures [39] and up to 0.25 eV at room temperature [40] 
opened in the bilayer graphene, the size of the gap being 
tunable by the electric field. 



On the other hand, according to density functional 
theory (DFT) calculations, cobalt adatoms on graphene 
show the highest magnetic anisotropy energy relative to 
other transition metals adsorbed on graphene [41, 42]. 
Thus, the Co-graphene system has potential applications 
in magnetic data storage. Recently, scanning tunnel- 
ing microscopy (STM) and spectroscopy measurements 
demonstrated that individual Co atoms deposited onto 
back-gated graphene devices can be controllably ionized 
by application of a gate voltage or STM tip bias volt- 
age [43]. Several different features were observed in 
the gate-dependence of the measured differential conduc- 
tance. Inelastic electron tunneling and hybrid electronic- 
vibrational states were suggested as possible explanations 
for some of these features. Therefore, the effects of in- 
elastic scattering of electrons should be included in any 
theoretical study whose objective is to reproduce these 
features of the calculated inelastic tunneling spectra [44- 
47]. A different feature is related to a sharp peak in the 
spectrum which moves in the opposite direction in en- 
ergy than the graphene Dirac point energy in response 
to changes in the gate voltage. 

In a previous theoretical study, using DFT electronic 
structure calculations in a generalized gradient approx- 
imation for periodic systems, Jacob and Kotliar [12] 
showed that a single Co atom most likely adsorbs at 
the hollow site of a clean graphene sheet. They found 
that the coupling of the Co 2>d levels to the graphene 
substrate and the dynamic correlations arc strongly de- 
pendent on the orbital symmetry, temperature, distance 
of the Co atom from the graphene sheet, and gate volt- 
age. More recently, Chan et al. [13] investigated the 
gated Co-graphene system in a superccU structure with 
periodic boundary conditions using the local density ap- 
proximation (LDA)-I-C/ method, with several values of 
the U parameter. The projected density of states was 
calculated for different U values and back gating treated 
qualitatively by LDA-f- U calculations for different inte- 
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ger charge states of the sample per unit cell. An ion- 
ization effect for the Co adatom was found as the charge 
state of the whole system was varied. No calculated STM 
current-voltage characteristics were reported, so the cal- 
culated local density of states (LDOS) could not be com- 
pared with STM data. Hence, the true LDOS of the Co- 
graphcnc system remains uncertain. Note that experi- 
mental and theoretical studies of the interaction between 
Co and graphene have shown charge-transfer to occur 
between the graphene and Co, resulting in the graphene 
being electron doped [48-50]. 

The aim of this paper is to reproduce some resonance 
features seen in the STM experiment by means of DFT 
calculations of the relaxed geometry and the extended 
Hiickel model of quantum chemistry [51-53] for the elec- 
tronic structure of Co adatom on graphene with DFT- 
based spin splitting parameters. To this end, we study 
resonant scattering of spin-up and spin-down electrons 
by the single Co atom adsorbed on gated graphene for 
electron energies close to the Dirac point. The electronic 
density of states for the Co-graphene system is calculated 
and the dependence of the resonant states on the gate 
voltage is discussed. In addition, we investigate the dif- 
ferential conductance of electron tunneling through the 
Co adatom as a function of sample bias and gate voltage 
which has not been addressed in the previous theories. 
Our transport calculations show that the experimentally 
observed behavior of the previously unidentified S-pcak 
in the STM differential conductance spectrum [43] is con- 
sistent with this feature being due to tunneling between 
the STM tip and a cobalt-induced Dirac point resonant 
state of the graphene, via a Co 3c? orbital. In this study, 
we do not consider the Kondo effect of Co on graphene 
[9, 10, 12], because this is a many-body effect and outside 
of the scope of the present study. 

The paper is organized as follows. In Sec. II, 
we present the model Hamiltonian of the infinite 2D 
graphene sheet with a single cobalt adatom. We show 
that the tight-binding Hamiltonian for the system can 
be transformed into an effective graphene Hamiltonian 
that includes the effect of the Co adatom. In Sec. HI, 
we derive an exact relation for the graphene Green's 
function matrix elements in the presence of the adatom 
using T-matrix theory. The results for the spin- and 
gate-dependent T-matrix and the electronic density of 
states will be presented in Sec. IV, where the adsorbatc- 
induced resonance states are discussed. The effect of 
local rehybridization of the graphene from the sp^ to 
sp^ electronic structure on the Dirac point resonances, 
that occurs when the Co atom bonds covalently to the 
graphene, is also examined. In Sec. V, we present a 
formalism suitable for studies of the gate-induced ion- 
ization of a single Co adatom on graphene, seen in the 
experiment. We calculate the gate dependence of the dif- 
ferential conductance by means of the Landauer-Biittiker 
formalism based on the Green's function technique, and 
the extended Hiickel theory and DFT-based spin split- 
ting parameters and compare the results with the STM 



measurements of Brar et aZ.[43]. Finally, in Sec. VI we 
conclude this work with a general discussion of the re- 
sults. 

II. THE MODEL 

Our calculations are performed using ab initio geom- 
etry relaxation based on density functional theory for 
a Co adsorbate on the honeycomb graphene lattice us- 
ing the Gaussian 09 software package together with the 
HSEhlPBE hybrid function and the 6-311G(d) Gaussian 
basis set [51-53]. Such a relaxed geometry is expected 
to be accurate because the density functional theory on 
which this relaxation is based, has already been well opti- 
mized for carrying out accurate ground state total energy 
calculations. We found that the lowest-energy binding 
site is the H-site in which the adsorbate is positioned at 
height ft-o = 1-68 A above the middle of a hexagon of the 
graphene, as shown in Fig. 1. The reason is that, the 
H-site has the largest number of neighboring C atoms, 
lowering the adsorption energy. The structure studied 
was a graphene disk consisting of 54 carbon atoms pas- 
sivated at the edges with 18 hydrogen atoms and the Co 
atom being bonded to the graphene at the center of the 
disk. All the atoms were allowed to relax freely. Al- 
though a finite size cluster of carbon atoms was used in 
the calculations, the bond lengths between the Co atom 
and its neighboring C atoms were well converged which 
indicates that 54 C atoms is large enough to model the 
graphene plane. The final relaxed structure obtained in 
this way is shown in Fig. 1. From the figure, it is clear 
that, there is no distortion out of plane of the graphene 
structure due to the Co adsorbate. This is in agreement 
with the previous studies of Co-graphene systems in pe- 
riodic supercell structures [7, 13]. In addition, we found 
in-plane distortion of 0.01 A in bond length between car- 
bon atoms of central hexagon relative to their nearest 
neighbor carbon atoms in the graphene plane. The ex- 
istence of an in-plane distortion equal to 0.023 A in the 
periodic Co-graphene supercell has also been reported by 
Liu et al. [16]. 

We describe the present system by a tight-binding 
model Hamiltonian derived from extended Hiickel theory 
in a basis of extended molecular orbitals (EMOs) which 
are linear combinations of the atomic valence orbitals 
of the adatom and the 2s, 2^3,, and 2py orbitals of the 
graphene carbon atoms. We refer the reader to Ref. [29] 
for a detailed explanation regarding tight-binding models 
for adsorbates on graphene derived from extended Hiickel 
theory and to Refs. 54-56 for further discussions of the 
application of extended Hiickel theory to graphene-bascd 
systems and carbon nanotubes. Accordingly, the EMOs 
together with the 2pz orbitals of the graphene carbon 
atoms form the basis set for our total Hamiltonian of the 
Co-graphene system that can be written in the form 

) (1) 
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FIG. 1. (Color online) Side view of the relaxed geometry for 
a single adsorbed Co atom on graphene. The optimal height 
ho = 1.68 A is the perpendicular distance between the Co 
adsorbate and graphene plane, while /i = 4 A is the distance 
between the adsorbate and STM tip which is directly above 
the Co atom. The inset shows the top view of the Co-graphene 
system. The hydrogen atoms at the edges of graphene are not 
shown here. 

where the first term corresponds to the Hamiltonian of 
clean graphene in the basis of 2pz carbon orbitals, a„cr is 
the annihilation operator for an electron with spin a in 
the 2pz orbital (f>na of carbon atom n, d\.^ is the creation 
operator for an electron in an EMO ■i/'ao- of the Co adatom 
and eccr is the corresponding energy eigenvalue, ^an.a 
is the spin-dependent matrix element of the extended 
Hiickcl Hamiltonian between the 2pz orbital of carbon 
atom n and the EMO a. We note that, since the basis 
set used in the extended Hiickel theory is nonorthogonal, 
it is possible for the overlap San,a = {''i'a(j\4>na) between 
the 2pz orbital 0„cr of carbon atom n and EMO ■0ao- to be 
non-zero. Because of this, we replace 7an,(T in the Hamil- 
tonian by ^an,a — ^San,a where e is the electron energy 
[57, 58]. Henceforth we will omit the spin index a in the 
overlap parameter and show it as Sq,„ for simplicity. 

The unperturbed Hamiltonian of two dimensional 
graphene, Hq, which can be modeled on a honeycomb 
lattice of two nonequivalent sites (that is, A and B) per 
unit cell is expressed as [6] 

HQ^^tY,[a\a,+\,.c.] (2) 

{hi) 

and describes the motion of carriers between 2pz carbon 
orbitals of graphene with a nearest-neighbor hopping pa- 
rameter t. Wc neglect the small changes in t due to the 
above mentioned in-plane distortion in the graphene ge- 
ometry, induced by the Co atom. We note that although 
Eq. (2) does not include the overlaps between the nearest 
neighbor 2pz orbitals of the carbon atoms explicitly, it is 
a widely used form of the Hamiltonian that has been opti- 
mized to provide an accurate description of the graphene 
TT-band states around the Dirac point. [6] We adopt it in 
the present work for that reason. Furthermore, since it 



is assumed that in clean graphene the spectrum is spin- 
degenerate, for simplicity we have suppressed the spin 
indices in Eq. (2). Here a\ is the creation operator 
for an electron in 2pz carbon orbital i. The Hamilto- 
nian for graphene is diagonal when written in the basis 
|k,c(w)) defined as |k,c(w)) = (|k, A)±e-*'*('')|k,B))/^/2 
where A and B stand for A and B sublattices and 
1'"'^^^)^ = 7^ ^^MB, Here, the 

c {v) sign denotes the conduction (valence) band, |R^(s)) 
represents an atomic orbital centered at a lattice site 
A{B) and Na{b) is the total number of A(i3)-lattice 
points. The phase (/)(k) is the polar angle for mo- 
mentum k and is defined as (/)(k) = Im(ln/(k)) where 
/(k) = e**^-"/^ 2e-*'''-"/2^cos(fc.ya/2). The corre- 
sponding eigenvalues in terms of /(k) can be expressed 
as£,(„) -±tv^7(k)P [5, 6]. 

The carbon 2s, 2px, and 2py orbitals in graphene com- 
bine to form in-plane cr (bonding) and a* (anti-bonding) 
orbitals. In this study the 2s, 2px, and 2py orbitals of the 
graphene were shifted down in energy relative to their en- 
ergies in extended Hiickcl model so as to yield an energy 
gap between the graphene a and tt bands in agreement 
with the results of the ab initio band structure calcula- 
tions for graphene [62]. 

The electron eigenfunctions of the total Hamiltonian 
can be written in terms of 2pz carbon orbitals and the 
adsorbed atom EMOs. In this regard, it has been shown 
[21. 29] that there is a relationship between 2pz orbitals 
of graphene and the adsorbate orbitals, so that, the 
total Hamiltonian Eq. (1) can be replaced by an ef- 
fective Hamiltonian which consists of the graphene 2pz 
Hamiltonian with an energy-dependent potential on each 
carbon atom to which the adsorbed atom binds. In 
fact, the energy-dependent effective potential acts like 
a self-energy that an adsorbate induces on electrons in 
graphene. For a single Co adatom that bonds to six 
neighboring graphene carbon atoms, the effective Hamil- 
tonian can be written as [29] 

iJeff =Ho+V (3) 
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^ = X! X! ^nm,Ga^naO'ma (4) 
n.m — 1 a 

where ana is the annihilation operator for an electron 
in the 2pz orbital of carbon atom n and Vnm.a = 
Y.a lan,al*a„i,a / " ^acr) displays a rcsouaut energy de- 
pendence. Here, the summation is over the extended 
molecular orbitals a of the Co adatom that bonds to 
carbon atoms m and n of the graphene. In this study, 
we consider 171 EMOs for the adsorbed Co atom on 
graphene which are linear combinations of the Co Sd, 
4s and Ap valence orbitals and the 2s, 2px, 2py valence 
orbitals of 54 carbon atoms shown in Fig. 1. The spin de- 
pendence of and 7cm, CT originates from the spin split- 
ting of the Co 4s and 3d orbitals. Due to the crystal field 
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splitting, the Co 3d orbitals split into three symmetry 
groups: the singly degenerate Ai group which only con- 
sists of 3d3z2_^2 orbital, a doubly degenerate Ei group 
consisting of the Sdzx and 3dzy orbitals, and a doubly de- 
generate E2 group derived from the 3dxy and 3dj;2_y2 or- 
bitals. The spin splitting values appear as extra energies 
that modify the Co 4s and 3d orbitals on-site energies 
obtained from the extended Hiickel theory parameters. 
The value of these splittings are calculated from each of 
two major peaks in spin density of states spectrum for 
Co from Ref. [13]. In addition, since in the parameters 
obtained from the extended Hiickel theory the effect of all 
carbon atoms (neighbors) is included in the calculations, 
the Dirac point of graphene will shift somewhat (0.072t) 
from zero energy to higher energy. To compensate for 
this effect, the cobalt energy levels are also shifted in 
energy by the same amount. 



III. T-MATRIX AND UNPERTURBED 
GREEN'S FUNCTION 

To study the strength of scattering associated with a 
single Co adatom on graphene we use the T-matrix for- 
malism which is defined as = G° + G'^TG'^ where 
G^^ = {e + irj — Hes)~^ is the effective Green's function 
for the graphene in the presence of the Co adatom, and 
G'^ ~ {e + irj — Hq)^^ is the unperturbed Green's func- 
tion (i.e. graphene Green's function in the absence of Co 
adsorbate) for 7r-band electrons. The T-matrix is related 
to G° and V hyT ^V{1- G^V)'^ where all the opera- 
tors are 6x6 matrices with the basis of the graphene pz 
orbitals of the carbon atoms to which the Co atom binds. 
To obtain the T-matrix components we need the matrix 
elements of Go- Since, e«(k) = — ec(k) the matrix ele- 
ments of {i\G^{e)\j) = G°j{e) between two arbitrary sites 
i and j of a graphene sheet can be written as 



(z|k,+)(k,+|j) , (z|k,-)(k,-|j) 



e -I- i?7 - ec(k) e -I- 177 -(- ec(k) 



(5) 



where the sum is over the Brillouin zone, e is the energy 
and ?7 is a positive infinitesimal. The results for G^j (e) 
depend on whether the sites i and j correspond to two 
equivalent or inequivalent sites. In the case of equivalent 
sites, the diagonal and off-diagonal matrix elements can 
be obtained from the following equation 



k 



(e + — ec(k) 



ik-(Ri-Rj) 



(6) 



where N 



N 



MB)- 



For the case of inequivalent sites. 



the off-diagonal matrix elements of G" corresponding to 
two inequivalent atoms (A and B in the same cells or in 
different cells) are written as: 



k 



-£c(k) 



(e + iiiy 



^i[k-(R.-R,)+0(k)] 



We note that for the cases of adsorbates such as H and 
OH that bond to only one carbon atom and also for O 
adatom which bonds to two neighboring carbon atoms, 
by expanding the band structure of graphene close to 
the Dirac points (at the K and K' points in the Brillouin 
zone), one can linearize the dispersion relation and obtain 
analytic expressions for the required matrix elements that 
is valid in the energy range \e\/t < 0.8 (sec Ref. [29]). For 
a general matrix element of the Green's function, how- 
ever, there is no such analytic expression, and we there- 
fore consider the full 7r-band structure of graphene and 
evaluate numerically the above Green's function matrix 
elements which are valid in all energy ranges. 

The locations of resonance energies are determined by 
the poles of T-matrix, which can be different for the spin- 
up and spin-down states. By taking matrix elements of 
the effective Green's function between the graphene 2pz 
orbitals of the six carbon atoms to which the Co atom 
binds, we obtain 



GlU^) = G%{e) + Gl{e)Tu,A^)Gl,,{^) (8) 
ii'=i 



where Tiii,„{e) are the spin- and energy-dependent T- 
matrix elements which describe resonant scattering of 
electrons in graphene due to the adsorbate levels renor- 
malized by hybridization with 2pz orbital of carbon 
atoms. The effect of this hybridization usually appears 
as a shift in the resonance levels towards the Dirac point 
of graphene [21, 29]. 



IV. SPIN-DEPENDENT DIRAC POINT 
RESONANCE OF CO ADATOM 

A. T-matrix feature 

To investigate the scattering effects associated with 
a single Co adatom on graphene, we study the square 
modulus of the T-matrix. The scattering will be a spin- 
dependent process due to the spin-splitting of Co 4s and 
3d orbitals, and the scattering at the energies e at which 
|Tcr(e) 1^ has maxima is of the resonant type. The energies 
of these resonant states can be found by setting the de- 
nominator of the square modulus of the T-matrix equal 
to zero, i.e. |det(l — G^V)\^ ~ 0. Based on the ^an,a 
and tan,a parameters obtained from the present model, 
the square modulus of the T-matrix for spin- up and spin- 
down electrons is shown vs. the electron energy e in Fig. 
2. The two studied models [solid and dashed curves) are 
explained in detail below. Note that, the effect of over- 
laps San between the EMOs ipa and the 2pz orbitals 0„ 
of carbon atoms n are included in all of the calculations. 

The solid curves correspond to a model in which all 
171 EMOs arc used in the calculations. In other words, 
the EMOs are linear combinations of the 4s, Ap^, 'ipy, 
(7) 4pz, 3da;2_j,2, 3dxy, 3(i322_^2, 3dzx, 3dzy atomic valence 
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FIG. 2. (Color online) Spin-dependence of square modulus of 
T-matrix vs. electron energy at the gate voltage for which 
Eg , the gate induced shift in energy of the graphene orbitals 
relative to those of the Co, is zero. Solid lines correspond 
to a case in which the EMOs are linear combinations of the 
atomic valence orbitals of the adsorbed Co atom and the 2s, 
2pj;, and 2py valence orbitals of each of the C atoms, while the 
dashed lines are related to a case in which only the valence 
orbitals of the adsorbate (no carbon orbitals) are included in 
the EMOs. The symbol t il) corresponds to spin-up (spin- 
down) electrons. The inset shows the effect of gate voltage on 
the square modulus of T-matrix when the 2s, 2px, and 2py 
orbitals of carbon atoms are included in the calculations. T^, 
e and Eg are in units of t—2.7 eV. 



orbitals of the adsorbed Co atom and the 2s, 2prr, and 
2py valence orbitals of all 54 carbon atoms. The square 
modulus of the T-matrix of the spin-down electrons dis- 
plays a resonant peak at the Dirac point of graphene, 
while for the spin-up electrons there is no such a reso- 
nance feature in the shown energy interval (|e|/t < 0.3). 
This is a measure of high spin polarization in the system. 

The dashed curves correspond to a simpler model of 
the EMOs in which no carbon orbitals are taken into ac- 
count in the EMO calculation and only 9 orbitals due to 
the Co atom arc included in the EMOs. We see that, the 
omission of the graphene carbon 2s, 2px and 2py valence 
orbitals in the calculations and of the coupling of adsor- 
bate to these orbitals (which are involved in the local 
rehybridization of the graphene from sp^ to sp'^ bond- 
ing near the adsorbed Co atom) shifts the position of 
the adsorbate-induced Dirac point resonance energy to 
eoR = — 0.053i and its intensity is about one order of 
magnitude weaker than that found in the more complete 



model. 

In earlier work [25, 29, 60, 61] it has been shown that 
the adsorption of H, F, OH, and O atoms on graphene 
distorts the graphene plane by shifting carbon atoms out 
of the plane. The omission of the graphene carbon 2s, 
2px, and 2py valence orbitals from EMOs was found to 
affect the adsorbate-induced resonances for these adsor- 
bates very strongly [29]. On the other hand, we have 
already explained that there is no out of plane structural 
distortion due to the Co adsorption on graphene, while 
if these valence orbitals are omitted from the EMOs ipa, 
the Co resonance state is still affected significantly. Since 
the inclusion of 2s, 2px, and 2py valence orbitals in the 
EMOs plays a dominant role for the Co adsorbate on 
graphene, and these orbitals play a crucial role in the 
rehybridization of graphene from sp^ to sp'^ bonding, it 
seems reasonable to infer that significant rehybridization 
from the sp^ to the sp^ electronic structure occurs for 
this system, despite the absence of structural out of plane 
distortion of the graphene. 

In the experiment of Brar et al. [43] both a back gate 
and the STM tip played the role of gates. Together they 
induced electrostatic potential differences between the 
graphene plane and the adsorbed Co atom that mani- 
fested themselves in the experimental STM data [43] . We 
model this gate-induced potential difference phenomeno- 
logically by shifting the carbon orbitals energies, derived 
from the extended Hiickel theory, and also the electronic 
structure of graphene 7r-band used in the unperturbed 
Green's function of the graphene plane, G", by a gate 
voltage-dependent amount Eq relative to the energies 
of the Co orbitals. Thus Eg acts as a proxy for the 
gate voltage in the present work. Accordingly, Eg = 
corresponds to the gate voltage at which the graphene 
Dirac point and the resonant state due to the adsorbate 
are located at zero energy e/t = 0.0. The gate volt- 
age at which Eg = is an important reference point 
for the gate-controlled ionization of the Co atom, espe- 
cially when we compare our differential conductance re- 
sults with the STM and spectroscopy measurements in 
the Sec. V. We emphasize that, in our calculations, the 
extended Hiickel model parameters describing both the 
on-site and intersite Hamiltonian matrix elements are re- 
calculated for each value of the gate voltage. The depen- 
dence of the intersite Hamiltonian matrix elements on 
the gate voltage [59] implies that hybridization between 
the orbitals on different atoms depends on the gate volt- 
age in the present theory. We note that previous work of 
Sofo et al. [60, 61] has indicated that doping can lead to 
a strong change in hybridization of the graphene. 

The application of a gate voltage in the experiment 
[43] caused the Co density of states features to move 
relative to the Fermi energy and hence, the Co levels 
could be emptied or filled with electron carriers. [43] The 
movement of a prominent peak (marked S in the STM 
conductance spectra shown in Figs. 2a and 2b in Ref. 
[43]) relative to the Dirac point with changing gate volt- 
age can be interpreted as further evidence of a Co en- 
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TABLE I. Minimal set of effective tight-binding parameters eaa, Tan.o- and s^n in units oit = 2.7 eV for adsorbed Co atom on 
graphene at zero gate voltage. The spin-dependent EMO energies e^o- are measured from the Dirac point energy of graphene. 
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-0.855 
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-0.044 
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-0.038 
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-0.001 
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-0.068 
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A OOO 

U.228 
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ergy level moving relative to the graphene tt band as the 
gate voltage was varied. The strength of the conductance 
peak S was observed to increase greatly in the experi- 
ment (Fig. 2a in Ref. [43]) as the peak approached the 
graphene Dirac point (Fig. 2b in Ref. [43]). As will be 
discussed in Section V, we interpret this behavior as be- 
ing due to a Co resonant level approaching the graphene 
Dirac point as Eq approaches zero. As this happens the 
cobalt-induced resonance feature for spin-down graphene 
electrons (associated with the peak in the square mod- 
ulus of T-matrix in Fig. 2) also approaches at the Dirac 
point as Eq strengthens. However with increasing \Eq\, 
the probability for overlap between the resonant cobalt 
level for spin-down electrons and the graphene states in 
the vicinity of the Dirac point decreases and the sharp 
conductance peak is suppressed in a similar way to the 
rather broad peak (marked S) in Fig. 2 in Ref. [43] for 
the higher absolute values of the gate voltage. 

The inset of Fig. 2 shows the calculated Eq depen- 
dence of the square modulus of the T-matrix and indi- 
cates that the electron energy ^ at which the resonance 
occurs and the strength of electron scattering depends 
on the value of gate voltage. We have mentioned above 
that a sharp resonance peak centered at cdr = 0.0 1 
appears at the gate voltage for which Eq = due to 
spin-down electrons. However, with increasing Eq the 
strength of resonance peak in the inset of Fig. 2 decreases 
and the peak position shifts towards the higher energies; 
ez3i?, = 0.045 1, 0.078 1, 0.11 1 at Eg = 0.2 0.4 0.6 < re- 
spectively. In the case of decreasing Eq the peak position 
cdr^ shifts to lower energies. For spin-up electrons the 
variation of the ]T£r(e)l^ function vs. energy is approx- 
imately the same in the all gate voltages. The reason 



of such a behavior can be understood in terms of LDOS 
of the Co adsorbate and graphene sheet that will be ex- 
plained below. 

The EMO energies Eqct, hopping jan.a and overlap San 
parameters for a minimal tight-binding model Hamilto- 
nian that describes the Dirac point resonant state in- 
volved in the spin-dependent square modulus of the T- 
matrix and graphene LDOS at zero gate voltage are pre- 
sented in Table L 



B. Local Density of States 

The calculated spin-dependent LDOS for graphene in 
the presence of an adsorbed Co atom, and for the Co 
4s, Ai, El, and E2 orbitals at different gate voltages 
are shown in Fig. 3. For greater clarity, the electronic 
states within an energy window close to the Dirac point of 
graphene are also shown separately [see Fig. 3(a2)-(d2)]. 
In the case of the graphene LDOS, we used the same 
extended Hiickel-theory-based model in the calculations 
that was used to calculate the square modulus of the T- 
matrix. Thus, the model includes the effects of the local 
rehybridization of the graphene plane due to the pres- 
ence of the Co adatom. The LDOS in this case is defined 
as A,a(e) = -lni[{i\Gf{€)\i)]/TT, where Gf{e) is given 
by Eq.(8) and \i) represents a 2pz orbital of a carbon 
atom. The solid curves in Fig. 3 correspond to the 2pz 
orbital of a graphene carbon atom to which the Co atom 
bonds. On the other hand, for the LDOS of the Co adsor- 
bate, we have used the Hamiltonian = + S^(e) 
for an electron with spin a (=t or l) where is the 
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FIG. 3. (Color online) Calculated LDOS per site vs. elec- 
tron energy e for a single Co atom adsorbed on grapliene 
at various gate voltages: (al), (a2) Eg = 0.0 i; (bl), (b2) 
Eg = 0.2 1, (cl), (c2) Eg = 0.4 i; (dl), (d2) Eg = 0.6 1. The 
symbol t (4-) corresponds to spin-up (spin-down) electrons. 
The LDOS associated with the 2pz orbital of a graphene car- 
bon atom to which the Co atom bonds (solid black), Co 4s 
(dash-dotted blue), Co 3d Ai (dash-dot-dotted dark yellow), 
Co 3d -El (dashed red), and Co 3d E2 (dotted green) states are 
shown. In Co LDOS, we multiplied by the scale factor 1/50. 
The EMOs included in graphene LDOS calculation are linear 
combinations of the atomic valence orbitals of the Co atom 
and the 2s, 2px , and 2py valence orbitals of C atoms. For Co 
LDOS the effect of carbon 2p~ orbitals is included as a self- 
energy in the calculations. The red and blue arrows indicate 
the position of the graphene Dirac point and the graphene 
resonant state, respectively. 



Hamiltonian of a free Co atom whose matrix elements are 
given by the extended Hiickel theory, and the self-energy 
S^(e) — Tco-c,<TG°(e)rQQ_Q ^ contains the coupling ef- 
fects of the adsorbate to the graphene carbon atoms [63] . 
Here, tco-Co- is a hopping matrix between the cobalt 
atom and the 2pz orbitals of the six carbon atoms of 
the graphene whose matrix elements are also given by 
the extended Hiickel model parameters and ^^(e) is the 



Green's function of the clean graphene which is given by 
Eqs. (6) and (7). In the above self-energy, the effect of 
2s, 2px, and 2py valence orbitals of C atoms could also 
be considered in the calculation. In such a case, how- 
ever, the TT-band Green's function, G°(e), is no longer 
sufficient and the graphene cr-band Green's function is 
required. Since the present self-energy is still able to 
reproduce the expected resonance features, we consider 
just the hopping parameters between the cobalt orbitals 
and the graphene 2pz orbitals in the Co LDOS calcula- 
tion. Note that the overlaps sco-c between the Co atom 
and the 2pz orbitals of the carbon atoms to which the 
adsorbate bonds are included in the calculations. Ac- 
cordingly, the LDOS for the adsorbate is determined by 
Di^ai^) = — Im[(£|G'(T(e)|.^)]/7r, where \£) represents the 
4s, Al, El, or E2 orbital, and Go-(e) = (e + ir] — Ha)~^ 
is the Green's function of the Co atom adsorbed on the 
graphene plane. 

From Fig. 3 it is clear that the spin-down graphene 
states show a resonance at the energies close to the Dirac 
point that is due to coupling between the graphene and 
the Co El orbital. This resonance is exactly positioned at 
the Dirac point energy at the gate voltage for which Eq, 
the gate induced shift in energy of the graphene orbitals 
relative to those of the Co, is zero. For this gate voltage 
the spin-down Co Ei state coincides with the resonant 
graphene state that appears due to the presence of the 
Co adatom. Furthermore, because of the spin-splitting 
of the Co 3d orbitals, the peak energy of the spin- up Co 
El state is somewhat far (i.e., at approximately — 0.55t 
for Eq = 0) from the Dirac point energy and accord- 
ingly, only a relatively broad peak due to the spin-up 
Co El state appears in the spectrum. The Co Ai and 
E2 orbitals are completely occupied for both spin states 
at the gate voltage for which Eq — 0. Note that, the 
Fermi energy Ep is set at the Dirac point of graphene, 
i.e. Ep = 0.0 1. In this study, only the effect of posi- 
tive Eq is investigated. The reason for this will be ex- 
plained in the next section. In both the graphene and 
Co LDOSs there is an energy shift with changing gate 
voltage. Although, as we mentioned above, the energy 
level shift due to the gate voltage is only applied on the 
graphene plane, the Co orbitals are also affected by this 
effect because, due to the non-zero overlap between the 
Co and graphene atomic orbitals, the effective hopping 
energies between the carbon and the Co orbitals depend 
on the on-site energies of the carbon orbitals [57, 58] . In 
the energy window {\e\/t < 0.3) that is shown in Fig. 3, 
the gate voltage decreases the occupation of both spin- up 
and spin-down graphene states by shifting the electronic 
states above the Fermi level Ep, and hence the occu- 
pation of spin-down Ei orbital is also decreased. This 
result can be clearly seen in Fig. 3 where Eq increases 
from O.Ot to 0.2t to 0.4t and 0.6t and the Dirac point 
of the graphene (indicated by red arrows) moves from 
zero energy to the values of 0.2t, O.At and 0.6t , respec- 
tively. Therefore, the Co orbitals which are coupled to 
the graphene shift somewhat to the higher energies and 
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the Co El orbital is gradually emptied of electrons. This 
result for Co 3d orbitals is in agreement with the recent 
DFT calculations for gated Co atom on graphene in su- 
percells with periodic boundary conditions [12, 13, 15]. 
Furthermore, by comparing Figs. 3(a)-3(d) it is clear 
that, with increasing Eg, the amplitudes of the Co Ei 
and graphene density of states peaks decrease, the en- 
ergy difference between the Dirac point and the graphene 
resonance state position increases and both the Co Ei 
level and the graphene resonance state induced by Co 
adatom are broadened. 

This indicates that the Co Ei orbital is coupled 
strongly to the 7r-orbitals of graphene and a gate volt- 
age can ionize the Co adatom through this orbital. This 
is in agreement with the work of Liu et al. [16] in which 
they reported that the hybridization between Ei orbital 
of the 3d-transition-metal adatoms on graphene and 
orbitals of the graphene carbon atoms is strong and re- 
sponsible for strong covalent bonding with graphene. By 
contrast, both the Ai and E2 orbitals are strongly lo- 
calized and far from the Dirac point, however, the Ai 
orbital hybridizes more weakly than the E2 orbital with 
the graphene states. Note that, comparison of the work 
of Liu et al. [16] with the present results shows some- 
what different coupling strengths of the Co E2 orbital 
with 2pz orbitals of the carbon atoms which may be at- 
tributed to the periodic nature of the structures (4x4 
Co-graphene supercell) studied in their calculations. It 
may also be relevant that density functional theory-based 
electronic structure calculations are known to have fun- 
damental limitations that can strongly affect the results 
that they yield for the coupling between adsorbed atoms 
and graphene; see Section VII of Ref. 29 for a dis- 
cussion and specific examples. In addition, the spin-up 
and spin-down Co 4s orbitals are far from the Ep; thus 
the occupation of these orbitals is zero at low gate volt- 
ages, although these orbitals hybridize strongly with the 
graphene states. This result is different than the results 
for periodic structures [7, 9, 13] in which the Co 4s orbital 
is close to the Dirac point energy and the effect of gating 
can ionize the Co atom by occupying or unoccupying the 
4s orbitals. In addition, in the periodic structures, the 
Dirac point resonances are broadened due to the presence 
of multiple adsorbed atoms on the graphene [7, 9, 13]. 
Therefore, it is evident that to better understand the fea- 
tures of adsorbate- induced Dirac point resonances and to 
verify the veracity of the theoretical predictions, compre- 
hensive experimental investigations would be very impor- 
tant. 



V. SIMULATION OF STM MEASUREMENTS 
A. Model 

In this section, assuming coherent transport, we 
describe how the electronic differential conductance 
(dl/dV) through a Co atom adsorbed on a graphene 



plane can be calculated using the tight-binding Hamil- 
tonians described in the preceding sections and the 
Landauer-Biittiker formalism [63]. The aim of this is 
to make a comparison with the STM differential con- 
ductance measurements in which the STM tip was held 
directly above the individual Co adatoms deposited onto 
a back-gated graphene device at 4.2 K [43]. To simu- 
late the experimental condition, we model the STM tip 
(shown in Fig. 1) by a one-dimensional chain of Pt atoms 
with nearest neighbor hopping integral t. The Fermi en- 
ergy is fixed at the center of the Pt chain band structure. 
For Pt atoms we use 6s, 6px, Qpy, &Pz, 5dx2-y2, 5dxy, 
5^322 _^2, ddzx, ^dzy atomlc valence orbitals. Since each 
Pt atom in the tip is described by 9 valence orbitals, 
then we model the tip by nine semi-infinite leads of or- 
thogonal atomic orbitals, one orbital per lead site, with 
spacing (periodicity) a between orbital sites. Each of the 
one-dimensional leads is decoupled from the others, thus, 
there is no hopping between different lead orbitals. Ac- 
cording to the above description, the Green's function of 
the Co atom coupled to the graphene carbon atoms and 
the STM system can be written as 

GAe,V) - [e+tv-H^°~^^{e)~^l™{e+eVr' , (9) 

where S^™(e) = Tco-STM,(T5STM(e)Tco_sTM,CT is a self- 
energy for an electron with spin a and describes the cou- 
pling of the adsorbate to the STM tip, Tco-STM,cr is the 
hopping matrix between STM and the Co adatom whose 
matrix elements are given by the extended Hiickel the- 
ory parameters, and gsTM = — (l/t)e'''°I is the surface 
Green's function matrix for the semi-infinite leads where 
k is the wave number in the leads and I is a 9 x 9 unit 
matrix [63, 64]. The overlaps sco-stm between the Co 
and Pt orbitals are included in the calculations. Now, 
if wc make use of the noncquilibrium Green's function 
technique to obtain the spin current la- for a sample bias 
V between graphene and the tip [63, 64], then the spin- 
dependent conductance dl^ / dV for an electron with spin 
a is given by 

GAV) liyA^,V) (10) 

where /sTM(e) is the tip Fermi distribution function, 
T,ie,V) = TT[T^ie)Gaie,V)rl™ie + eV)GUe,V)] is 
the Landauer transmission function that can be calcu- 
lated from the device Green's function, Eq. (9), and the 
broadening functions which can be expressed for elec- 
trode a (= C or STM) as F^ = -2Im(E^). In Eq. 10, 
the voltage dependence of transmission function has been 
ignored, because the transmission values are not sensitive 
to low value bias voltages. 

B. Results 

The differential conductance is the most direct way to 
measure the energy-dependent LDOS of a system. In the 
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present system, the calculated differential conductance, 
Q{V) = G^{V) + Gi{V) vs. sample bias V for different 
applied gate voltages is shown in Fig. 4. Note that, the 
sample bias, V, means the voltage of the graphene with 
respect to the tip. In the present study, the graphene 
remains at the zero potential and the tip is at negative 
voltages, similarly to the experiment. Hence, when a 
negative voltage is applied to the tip, the sample bias is 
positive {V > 0) [43]. Furthermore, the gate-dependent 
G{V) spectra in the experiment, were measured at nega- 
tive gate voltages. The measurements showed that, when 
the gate voltage moves to more negative values, the Dirac 
point position shifts to positive energies [see Fig. 2b in 
Ref. [43]]. This behavior is clearly equivalent to applica- 
tion of positive Eg values as we have done in Section IV 
and can be seen in Fig. 3. 

The calculated differential conductance G(V) shows a 
sharp peak at Eq = that comes from the Co spin- 
down resonant state. The height of this peak decreases 
while the spectrum is broadened as the gate voltage is 
varied and \Eg\ increases. This effect is due to the na- 
ture of resonant states induced in the graphene near the 
Dirac point by the presence of an adsorbate; as can be 
seen in Fig. 3 these resonant states become narrower and 
stronger as the associated energy level of the Co adsor- 
bate approaches the Dirac point energy, i.e.. Eg — > 0. 
When \Eg\ increases, the cobalt orbital goes out of res- 
onance with the graphene density of states peak in the 
vicinity of the Dirac point and the sharp conductance 
peak is suppressed. 

In the experiment of Brar et al. [43], it has been 
demonstrated that the difference between the location 
of Dirac point voltage (marked Vd in Fig. 2b of Ref. 
[43]), associated with the Dirac point energy [65], and 
the peak in the differential conductance signal (marked 
S in Fig. 2a and 2b of Ref. 43) increases as the mag- 
nitude of experimental gate voltage is increased. This 
enhancement in the difference between the Dirac point 
energy and the energy of the S state with increasing the 
gate voltage is in good agreement with the result of Fig. 3 
which shows that, although with increasing Eg both the 
cobalt-induced resonance feature (blue arrows) and the 
Dirac point energy (red arrow) in the LDOS spectrum 
of Fig. 3 move in the same direction, the Dirac point 
location moves more rapidly and the difference in energy 
between both locations (red and blue arrows) increases 
with Eg- 

The experimental data [43] exhibits a sharp resonant 
peak S in the G{V) spectrum at low gate voltages in Fig. 
2a of Ref. 43. With increasing absolute value of the gate 
voltage, the height of the peak decreases and the peak 
moves to lower sample bias. The sharp peak seen for 
Eg = in our model (Fig. 4), is similar to the observed 
resonant peak S in Fig. 2a of Ref. 43 at low gate voltages 
and evolves with gate voltage in a qualitatively similar 
way as well. 

The peaks in the GiV) spectra in Fig. 4 occur at the 
sample bias voltages at which the STM tip Fermi en- 



ergy crosses the peaks of the partial densities of states 
that are due to the Co d Ei state and the associated Co- 
induced graphene resonant state that are seen in Fig. 3 
a2-d2. These density of states features in Fig. 3 a2-d2 
and the peaks in GiV) spectra in Fig. 4 are strongest and 
narrowest when the energy of the Co d Ei state coincides 
with the graphene Dirac point energy, which occurs for 
Eg = 0. The same is true of the Dirac point resonance 
peak of the T-matrix as can be seen in Fig. 2 and its inset. 
The large increase in the intensity of the experimentally 
observed resonant peak S in Fig. 2a of Ref. 43 (with 
decreasing absolute value of the gate voltage) also occurs 
as the S-peak energy approaches that of the graphene 
Dirac point energy (labelled Vd) in Fig. 2b of Ref. 43. 
We conclude that the large increase in the intensity of the 
resonant peak S in Fig. 2a of Ref. 43 with decreasing ab- 
solute value of the gate voltage is the experimental STM 
signature of the strengthening and narrowing of Dirac 
point resonance induced in the graphene by the Co (and 
of the associated Co level as well) as the Co level crosses 
the graphene Dirac point energy. Thus we propose that 
STM is able to detect the presence of adsorbate-induced 
Dirac point resonances in the graphene by monitoring 
the intensity of conductance peaks associated with or- 
bitals of the adsorbate as those adsorbate levels cross the 
graphene Dirac point with changing gate voltage. Sys- 
tematic experimental studies of adsorbate-induced Dirac 
point resonances in graphene based on this finding would 
be of considerable interest since experimental studies 
of these resonances have, to our knowledge, not been 
reported previously despite the considerable theoretical 
interest [7-29] in this topic. 

The gate dependence of G{V) also indicates that the 
value of Eg determines the ionization state of the Co 
atom, because by changing the value of Eg, the Co 3c? 
states can cross the Fermi energy and accordingly the 
charge occupation in these states changes. Therefore, 
the application of gate voltage enables tuning the charge- 
carrier density of graphene and hence the variation of Co 
adatom ionization state. 

In the present model, the distance between the tip and 
the Co atom is h—A A. Therefore, the overlap and hence 
the direct coupling between the Co atom and the tip 
is weak and that between the tip and the graphene is 
much weaker still. The reason is that, the sum of Ft 
and Co atomic radii is ~ 2.64 A which is considerably 
less than h which results in a small overlap between their 
atomic wave functions. In reality, the Co levels couple the 
STM tip to the graphene to open a conductive channel for 
electron transport. Since the coupling between the Ei- 
levels of adsorbate and the STM tip is weak the transport 
process for such levels is tunneling and the differential 
conductance spectra are quite similar to the Ei resonance 
state LDOS seen in Fig. 3. From the extended Hiickel 
theory hopping parameters between the Co atom and 
the Pt atom, we found that the s-levels of both atoms 
hybridize more strongly than that the other levels. The 
hopping energy between s-levels is roughly one order of 
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magnitude greater than the hopping energy between Ei- 
levels. Note that, in this work because of the position of 
the Co i?i-level relative to the graphene Dirac point and 
also to the other levels such as Co 4s-level, the Ei level 
plays the main role in the tunneling transport through 
the cobalt atom at low values of the gate voltage. 



Finally, it should be mentioned that we also exam- 
ined the electron transmission through the Co adatom 
in the graphene-Co-Pt structure using the Landauer the- 
ory and the Lippmann-Sehwinger equation [66] in which 
the EMOs included in the graphene were linear combi- 
nations of the atomic valence orbitals of the Co atom 
and the 2s, 2pj:, and 2py valence orbitals of graphene 
carbon atoms. The results, except with a small shift in 
the energy of the resonant states, were the same that we 
obtained based on the present method. This confirms 
that by the inclusion of just carbon 2pz orbitals in the 
graphene self-energy, we can obtain the same resonance 
feature in the conductance spectrum that we obtained 
based on the Lippmann-Sehwinger approach. 



VI. DISCUSSION 

In this study, we have presented a tight-binding theory 
of the Dirac point resonances due to a single Co adatom 
on gated graphene based on the 7r-band tight-binding 
description of the graphene plane with nearest neigh- 
bor hopping energy t and the extended Hiickel model 
of the adsorbate and the adsorbate-induced local sp'^ re- 
hybridization of the graphene. Using the generalized ef- 
fective Hamiltonian [29] for the relaxed geometry of Co- 
graphene system based on ab initio density functional 
theory, we found a strong spin-dependent scattering res- 
onance at the Dirac point of graphene and a gate de- 
pendence of this resonance. This resonance arises from a 
hybridization of Co atomic orbitals with graphene states 
and such resonances have been predicted for 3d transi- 
tion metals adsorbed on graphene [7, 9, 11-13]. The en- 
ergy position of the Co atomic levels with respect to the 
graphene Fermi energy determines the ionization state 
of the adsorbate. Because of spin-splitting of Co 4s and 
3d orbitals, this ionization is different for spin-up and 
spin-down states. The results showed that, the Co Ei 
state can be emptied with electrons as the state shifts 
above the Fermi energy through application of gate volt- 
age, in agreement with the recent experiment [43]. The 
change in the occupation of the Co 3d orbitals and hence 
in the ionization state of the adatom by application of 
gate voltage has also been predicted using DFT-based 
calculations for periodic systems. [12, 13]. 

Furthermore, we presented a formalism that we used to 
model the differential conductance measurements. Our 
calculation shows that, when the Co adatom states cross 
the Dirac point, a sharp feature is induced in the GiY) 
spectrum which is consistent with experiment and has 
not been addressed theoretically, previously. The height, 
broadening and the position of this feature is strongly 
dependent on gate voltage. The agreement between this 
theoretical dependence and the experimentally observed 
behavior of the S-peak in the STM differential conduc- 
tance spectrum [43] allows us to attribute the S-peak 
to tunneling between the STM tip and a cobalt-induced 
Dirac point resonant state of the graphene, via a Co 3d 
orbital. Systematic experimental studies of adsorbate- 
induced Dirac point resonances in graphene based on this 
finding would be of interest. 

Our results also show the Dirac point energy of 
graphene and the cobalt-induced resonance feature in the 
differential conductance spectrum to move in different 
ways when a gate voltage is applied. Therefore the ion- 
ization state of Co adatom can be tuned by application of 
a gate voltage which induces a change in the occupation 
of the Co 3d states. From the above findings regarding 
the T-matrix and the electronic density of states we can 
conclude that the Co adatom on graphene has a strongly 
spin-polarized LDOS especially at electron energies in the 
vicinity of the Dirac point of graphene. This suggests 
that Co adatoms on gated graphene may have possible 
applications in magnetic data storage and magnetic sen- 
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sors and other applications which may be discovered. and Compute Canada. 
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